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Abstract 

We discuss several general aspects of the free energy of the standard model at high 
temperatures. In particular the Clausius-Clapeyron equation is shown to yield a 
relation between the latent heat and the jump in the order parameter. The free 
energy is calculated as function of temperature in resummed perturbation theory to 
two-loop order. A new resummation procedure is proposed in which the symmetric 
phase and the Higgs phase are treated differently. A quantitative description of 
the phase transition is achieved for Higgs masses below ~ 70 GeV. The results 
are found to be in agreement with recent numerical simulations on large lattices. 
The phase transition provides no evidence for strong non-perturbative effects in 
the symmetric phase. 
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1 Introduction 



In the standard model of electroweak interactions all masses are generated by the Higgs 
mechanism. At high temperatures this implies, in analogy to superconductivity, a phase 
transition from a massive low-temperature phase to a massless high-temperature phase, 
where the electroweak symmetry is "restored" jjj]. This transition is of great cosmological 
importance because baryon-number violating processes fall out of thermal equilibrium 
below the critical temperature of the phase transition 0. As a consequence, the present 
value of the baryon asymmetry of the universe has finally been determined at the elec- 
troweak transition. 

In recent years quantitative studies of the electroweak phase transition have been 
carried out by means of resummed perturbation theory 0]-|]7| and lattice Monte Carlo 
simulations P|-||12||. As a first step towards a treatment of the full standard model, the 
pure SU(2) Higgs model is usually investigated, neglecting fermion effects and the mixing 
between photon and neutral vector boson. These can be treated in perturbation theory 
and do not affect the essential infrared problem. There is general agreement that the 
phase transition is first-order for Higgs masses m# below the vector boson mass mw 



At larger Higgs masses lattice results fI5| , general arguments |14] and non-perturbative 
solutions of gap equations [T5j suggest that the first-order transition changes to a smooth 
crossover. However, this conjecture still remains to be firmly established. 

The goal of the present paper is a quantitative study of the thermodynamics of the 
transition for Higgs masses below 80 GeV. Our analysis will be based on the gauge 
invariant "order parameter" $t$ |Tj| and the corresponding free energy |17 |. We shall 



extend our previous one- loop analysis [18| to two-loop order and discuss the connection 



with the conventional Landau gauge approach. The comparison between one-loop and 
two-loop results, and also between different resummation procedures, will allow us to 
estimate the uncertainty of perturbative predictions for thermodynamic observables. 

It is well known that the electroweak phase transition is influenced by non- 
perturbative effects whose size is governed by the confinement scale of the effective 
three-dimensional theory which describes the high-temperature limit of the SU(2) Higgs 
model. These effects are particularly relevant in the symmetric phase, and one may worry 
to what extent a purely perturbative analysis of the phase transition can yield sensible 
results at all. However, it is not yet known how important these non-perturbative effects 
are quantitatively for different Higgs masses. As we shall see, at two-loop order our re- 
sults depend only logarithmically on the infrared cutoff needed in the symmetric phase. 
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We shall then compare the perturbative results with recent non-perturbative results ob- 



tained by numerical Monte Carlo simulations on large lattices [11]. This will enable us 
to estimate the size of possible non-perturbative corrections. 

The paper is organized as follows. In sect. 2 we shall discuss several general aspects 
of the free energy and the various effective potentials used in connection with the elec- 
troweak phase transition. The description of the transition is based on a gauge invariantly 
coupled source term. The corresponding gauge invariant free energy can be obtained from 
the minima of the usual Landau gauge potential. 

In sect. 3 a useful relation is derived between the latent heat and the gauge invariant 
order parameter, which follows from the Clausius-Clapeyron equation by dimensional 
arguments. This relation is exactly satisfied in perturbation theory, and it holds for the 
lattice results within the estimated errors. 

Sect. 4 deals with the free energy at one-loop. A new resummation procedure is 
introduced which treats the Higgs phase and the symmetric phase differently, reproducing 



the numerical results of |L8| in a different manner. This also clarifies the connection 
between the gauge invariant and the usual Landau gauge approach. It is shown that the 
barrier of the gauge invariant effective potential is given by analytic continuations of the 
convex potential defined by a Legendre transformation. 

In Sect. 5 higher order corrections are obtained using the new resummation method 
of the previous section. This calculation is an extension of the gauge invariant approach 



of |18| to the two-loop level. While previous numerical results for phase transition param- 
eters are qualitatively reproduced, the convergence of perturbation theory is improved 
in the new approach. 

Sect. 6 contains a comparison with lattice results, showing quantitative agreement 
within the statistical and systematic errors. For this the most important zero-temperature 
renormalization effects have to be included at small Higgs masses. 

Conclusions are summarized in sect. 7, and the appendix contains the explicit for- 
mulae used in sect. 5. 

2 Free energy of the SU(2) Higgs model 

The action of the SU(2) Higgs model at finite temperature T reads 

S>[$, W}= [ dx Tr V-W^W^ + (D^)W^ + + 2A($ t $) i 

J (3 L ^ 
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with 

$ = ^(<7 + i7?-f), D^ = (d li -igW^, W^U-W^ (2) 

, dx = I dr [ d 3 x , 13 = 1= . (3) 
p Jo Jn T 

Here H^ is the vector field, a is the Higgs field, tt is the Goldstone field and f is the 
triplet of Pauli matrices. In general, we shall consider the limit of infinite spatial volume 
Q. For perturbative calculations gauge fixing and ghost terms have to be added to the 
action ([!]). 

The free energy density of the system, W(T, J), is given by the partition function, 
i.e., the trace of the density matrix, 



exp(-pnW(T, J)) = Trexp 



-P (h + J J Q d 3 x&® 



(4) 



where H is the Hamilton operator of the theory, and $ is the operator describing the 
Higgs field. We have added a source J, with <9 M J = 0, coupled to the spatial average of 
the gauge invariant composite operator <J>t$ (here and below the trace operator acting on 
$t$ is omitted for brevity). Similarly, one may define a free energy for spatially varying 
sources J(x). The partition function can be expressed as a euclidian functional integral 



19 . 

exp(-pSlW(T, J)) = J D&D&DWpexp (~ J p dx ( L + ^ f $)) , (5) 

where L is the euclidean lagrangian density, and the bosonic fields $ and satisfy 
periodic boundary conditions at r = and r = /3. Eq. @ is the starting point of 
perturbative as well as numerical evaluations of the free energy. 

Note, that the source J in eq. (H) couples to a gauge invariant composite field. Hence, 
the free energy W(T, J) is gauge independent. The spatially constant source J simply 
redefines the mass term in the action (JJ). This is in contrast to the usually considered 
generating function of connected Green functions at zero momentum, 

exp (-pnW(T,j;J)) = J D&D&DWpexp J p dx ( L + J$t$ + J a )) ■ ( 6 ) 

Here the source j couples to a gauge dependent quantity, the field a. Consequently, 
W(T,j; J) is gauge dependent and not a physical observable. For later use we have also 
kept the dependence on the source J. 

The generating function W(T,j; 0) can be made finite in the usual way by a multi- 
plicative renormalization of couplings and fields. This is not the case for the free energy 
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W(T,J) } since J couples to a composite field. It is known [20] that two more counter 
terms, linear and quadratic in J, are necessary in order to subtract the additional diver- 
gencies. Hence, the renormalized free energy W(T, J) contains two arbitrary constants 
in addition to the usual renormalized parameters at zero temperature. 

From the free energy W(T,J) a gauge invariant effective potential V(T,p) can be 
obtained as usual by means of a Legendre transformation, 

V(T,p) = W(T,J)-±pJ , (7) 



where 



\p = 1 / d*x < &(x)$(x) >= %-W{T, J) 
2 \l Jo. oJ 



is the spatial average of the thermal expectation value of which plays the role of 
an "order parameter" in the SU(2) Higgs model. By definition, the effective potential 
V(T, p) is convex. The ground state of the theory corresponds to a stationary point 
of the effective potential, where dV(T,p)/dp vanishes. In the case of a first-order phase 
transition two stationary points connected by a straight line (see fig. 3 in sect. 4) represent 
two coexisting phases. 

What is the effect of the ambiguity of the renormalized free energy W(T, J) on the 
effective potential? Consider two definitions of the free energy, related by 

W(T, J) = W(T, J) + bJ + cJ 2 . (9) 

The two corresponding effective potentials are V(T,p) and 

V{T ) p) = W{T ) J)- l -pJ , (10) 

with 

l -p = ^W{T,J) = l -p + b + 2cJ . (11) 
For stationary points po and po of the effective potentials, where 

j-V{T,p) \ p=p = , ^V(T,p) \- p=Po = , (12) 

one easily verifies 

Po = p + 2b, V(T,po) = V(T,po) . (13) 

This means that the free energy of the ground state is independent of the two parameters 
b and c, whereas the expectation value p is arbitrary. However, in the case of more than 
one stationary point, which is relevant for first-order phase transitions, the difference Ap 
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between two stationary points is independent of b and c. This difference is a physical ob- 
servable. Other properties of the gauge invariant potential V(T, p), such as the curvature 
at a local minimum, will in general depend on the parameters b and c. 

The functional integral for the free energy may be written as 

exp (-P0W(T, J)) = 2n 2 J ip 3 d<p DQD&DWpS i^p - J^dx a^j 5 3 jT dx tt^J 

exp (- J dx(L + J$t$)) . (14) 

Here the integrand of the ordinary integral over <p is the exponential of the well known 
constraint effective potential |21], |22fl , 

««p(-/MW(r, w j)) = /wMt^^i-lt.j^i-li,) 

exp (-J dx(L + J$ f $)) . (15) 
In the infinite volume limit the constraint effective potential U(T, p\ J) coincides with 



the effective potential V(T, <p; J). From eqs. (|T4D,(l5D one obtains for the free energy, 

exp (-/3SW(T, J)) = 2vr 2 J p 3 dip exp (-/3QV(T, <p; J)) . (16) 

In the infinite volume limit, this yields 

W(T,J) = V(T,p min (T,J),J) , (17) 

where p m i n (T, J) is the global minimum of the effective potential V(T, ip; J). For arbi- 
trary values of <p the potential V is gauge dependent. However, its value at the minimum 
is known to be gauge independent [p3fl , yielding a gauge independent free energy W(T, J). 

At the critical temperature T c of a first-order transition the order parameter p and 
the energy density E(T,0), 

E(T, J) = W(T, J) - T^W(T, J) , (18) 

are discontinuous. The jump in the energy density is the latent heat AQ. This discontinu- 
ity in p and E requires a free energy W(T, J) with the following properties: W(T, 0) must 
be continuous but not differentiable at T = T c , and the same must hold for W(T C , J) at 
J = 0. In the following sections we shall verify these features based on a perturbative 
evaluation of the free energy. 
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3 Clausius-Clapeyron equation 



For the first-order phase transition from liquid to vapour there exists a well known 
relation between the latent heat and the change of the molar volume, the Clausius- 
Clapeyron equation EM. In the electroweak phase transition the "order parameter" 
plays the role of the molar volume and a completely analogous relation can be derived. 

The electroweak plasma can exist in two phases, the massive low-temperature Higgs 
phase with free energy Wb(T, J) and the massless high-temperature symmetric phase 
with free energy W S (T, J). In the J — T-plane the boundary between the two phases is 
determined by the equilibrium condition 

W s (T 1 J(T)) = W b (T 1 J(T)) , (19) 

which implies 

Using the definitions for latent heat and jump in the order parameter, 

-^(W s -W b ) = ~AQ , ±(W s -W b ) = ±Ap, (21) 

one obtains 

AQ = \^T% . (22) 
This is the Clausius-Clapeyron equation of the electroweak phase transition. 

So far, we have only used the continuity of the free energy along the phase boundary 
in the case of a first-order transition. We can now employ the fact that the mass term 
fi + J is the only dimensionful parameter of the SU(2) Higgs model. This implies 

// + J(T) = C(g 2 , \)T 2 , (23) 

and therefore 

J T(T\ 

1\i . (24) 



T dJ(T) 



dT j=o 

To leading order in the couplings one has (cf. ([34])), C(g 2 , A) = — (j^g 2 + |A). Inserting 
eq. (0) into eq. ( ^2j) we finally obtainf] 

AQ = -~m 2 H Ap(l + (D(g 2 ,\)) . (25) 



The higher order corrections are due to the difference between the mass parameter 

and the physical Higgs mass. In the following sections this relation will provide a very 

useful check on our results. 



lr The same result has been derived in a recent paper by Farakos et al. |25| based on renormalization 
group equations for the effective three-dimensional theory. 
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4 Free energy in perturbation theory 



Near the ground state, J = 0, the free energy W (T, J) can be evaluated as power series 
in the couplings g and A by means of resummed perturbation theory^. Here thermal 
corrections are added to the tree-level masses of the scalar fields and the longitudinal 
component of the vector boson field, 

8Sp = pj d 3 x {^a 01 T 2 (a 2 + vr 2 ) + \a{T 2 wt) • (26) 

The sum of tree-level masses and thermal corrections then enters the boson propagators 
in loop diagrams, and 8S% = —SSp is treated as counter term. In eq. (|26|) the fields a 
and 7r do not depend on the imaginary time r. Following Arnold and Espinosa [|j], we 
only resum the static modes of scalar and vector fields. Clearly, this is sufficient to avoid 
infrared singular contributions from these fields since non-static modes have thermal 
masses 0(T). For the usual resummation with one-loop thermal counter terms ~ T 2 
the resummation of static modes is known to be equivalent to the resummation of all 
modes. However, in the following we shall employ counter terms which depend on the 
scalar background field. It turns out that in this case the resummation of static modes 
is preferred, since only for this resummation the loop expansion is an expansion in the 
couplings g and A. 

To leading order in the couplings, one obtains for the parameters in eq. fl2T)|) from 
one-loop self energy corrections ||, 

«oi = y^9 2 + ^A , ai= ^g 2 . (27) 

The masses of the boson propagators are obtained from eqs. ([!]) and (p6|) by shifting the 
Higgs field a by the average field <p (cf. (|14|),( |l5|) ). This yields m L , m T , m a and for 
longitudinal and transverse part of the vector field, the Higgs field and the Goldstone 
boson field, 

m 2 L = a 1 T 2 + 9 ^ 2 , ml = g l^ i (28) 

ml = a m T 2 + n + J + 3A<^ 2 , m 2 = a 01 T 2 + /x + J + \<p 2 . (29) 

The scalar masses agree with derivatives of the effective potential V, 

ml = ^V{T^-J) , m2 = ~y(T l¥ ,;J), (30) 



2 For a detailed discussion and references, see 
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up to terms of higher orders in the couplings g and A. 

It has been shown that the resummed loop-expansion for the free energy is a sys- 
tematic expansion in the couplings g and A || . So far, however, several aspects of this 
resummation have remained unsatisfactory. First, higher order corrections are very im- 
portant for the scalar masses^. In particular, in the Higgs phase the Goldstone boson 
mass m n , as given by eq. (HO), vanishes at the minimum. Hence, the expression to leading 
order in the couplings, given in eq. fl29f ), is cancelled by higher order corrections. In the 
following we shall therefore modify the loop expansion in the Higgs phase. We replace 
eq. (H) by 

ml = 2Xcp 2 , m*=0, (31) 

and we treat 

5S b p = (3 J d 3 x^ (a 01 T 2 + n + J + \<p 2 ) (a 2 + n 2 ) , (32) 

as counter term, where a and ir again represent static modes. Note, that in the sum 5Sp + 
8S\ the temperature dependent terms cancel. Hence, we perform no thermal resummation 
for scalar masses in the Higgs phase. 

In the symmetric phase, <p = 0, and at temperatures close to the barrier temperature 



Tb = y/i/ctoi, the one- loop finite-temperature scalar masses are small, and higher order 
corrections are important for the Higgs boson and the Goldstone boson mass. Here, we 
will replace eqs. ( |2"U|) by self-consistently determined scalar masses, which are defined by 

ml = ml = -^-V(T,0;J) . (33) 

if oip 

For given vector boson masses and mj, this is a gap equation for the scalar masses, 
which can be solved at each order of the loop expansion. 

The self-consistent scalar mass at one-loop is easily obtained from the one-loop ef- 
fective potential which reads, for arbitrary resummations of the static modes, 

V{T,ip;J) = ?lj + ha 01 T 2 + l x + J)^ + -\^ 
o / 4 

" tL ( 3m ^ + 6m ^ +m ° + 3m ") + ° ( ^ 4 ' a2) • (34) 

From eqs. (|28"D , ( |33"D and (|34| ) one obtains in the symmetric phase 

aoT 2 + /i + J + O {g 2 mT, Xm^j , (35) 



2 2 



3 We thank P. Arnold for emphasizing this problem. 
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where 

*° - t/ + \ x - wA° 3 ■ < 36 > 

The corresponding counter term for perturbation theory in the symmetric phase, to be 
inserted in one-loop graphs in a two-loop calculation for W(T, J), reads 

SS' p = -(3 J d 3 x (~a T 2 (a 2 + tt 2 ) + \a{T 2 wt) . (37) 

Having specified the vector boson and scalar masses in the Higgs phase and the 
symmetric phase, as given in eqs. (|28|), (|3T| ) and fl35|), the free energy can be calculated 
from the effective potential V, using the relation (pTT|) . The effective potential has two 
local minima, 

Wmin,l{T, J) = (f s = , (Pmin,2(T, J) = (fib > , (38) 

which correspond to the symmetric and the Higgs phase, respectively. From eq. (0) one 
obtains for the free energy in both cases 

W S (T,J) = ^J-^( ao T 2 + /i + Jf 2 , 

O 671 

W b (T,J) = V(T,if b (T,J);J) (39) 

T 2 1 1 
= yJ + -(«oiT 2 + /i + J)ip 2 b + -Xipt 

~ ( 6 (l^^ 2 ) 3/2+ ( 2A ^) 3/2+ 3 + I^^) 3 " 3 «? /2t3 )-( 4 °) 

For convenience, we have subtracted the terms independent of J from W s and W b . 
The free energy of the ground state is given by 

W(T,0) =min{H/ s (T,0),H/ b (T,0)} . (41) 

It is a concave function, shown in fig. 1 for some choice of g and A. The derivative of 
W(T, 0) has a jump at the critical temperature T c , which is characteristic for a first-order 
phase transition. The corresponding latent heat is given by 

AQ = T A (W b (T, 0) - W S (T, 0)) \ T=T = (y 2 h + ^ (a T 2 + /i) 1/2 ) , (42) 

which follows immediately from the Clausius-Clapeyron equation or, with some more 
work, from eqs. (|39D , (|40|) . We can also consider W(T C , J), i.e., the dependence of the free 
energy on the external source at the critical temperature T c . This function, shown in fig. 
2, is also concave and similar to the function plotted in fig. 1. In this plot the huge linear 
term T 2 J/6, which has only the effect to shift the field square expectation value by T 2 /3, 
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has been discarded in both phases. As discussed in sect. 2, the jump in the derivative 
at J = yields the jump in the "order parameter" p at the first-order phase transition. 



One easily verifies that the Clausius-Clapeyron equation (25) is satisfied. 

From the free energy W(T, J) one can obtain the gauge invariant effective potential 
V(T,p) by means of a Legendre transformation. Since the derivative of W(T, J) is not 
continuous everywhere, one has to use the definition^], 

V(T,p)=su P {W(T,J)-l-pJ} . (43) 
J 2 

This yields the convex, non-analytic function, plotted in fig. 3 as full line. One may 
also compute the ordinary Legendre transform V s (T,p) and Vb(T,p) of W S (T,J) and 
Wb(T, J), respectively. Neglecting constant terms this yields 

V s (T,p) = IW + rip'-^p' 3 , 

V b (T,p) = V(T,Jp~',0), (44) 



where we have used 



1,1/ T 2 \ dW b (T,J) T 2 1 2 dV 



d(p b 1 2 



dJ 2 

ft 



-M • (45) 



V s and V b are also shown in fig. 3. In the region outside of the two local minima, V s and 
V b , respectively, coincide with the convex effective potential V(T,p). Between the two 
local minima, V s and V b represent two analytical continuations of V(T,p), which meet 
at the "matching point" pu = T 2 /3. At this point, marked by a cross in the plot, the 
first derivatives of both curves coincide. 

The non-convex "effective potential" obtained by combining V s and Vb on both sides 
of the "matching point" is almost identical with the gauge invariant "effective potential" 
obtained in fT8| . If the root appearing in Vb(T, p) is expanded to order p 2 /T 4 , which 



corresponds to the reduction to the three-dimensional theory considered in [JTSJ] , both 
"effective potentials" are identical. The generation of a barrier between two local minima 
as analytic continuation from a convex effective potential is reminiscent of the treatment 
of first-order phase transitions in condensed matter physics |[26|| . However, the precise 
physical meaning of the resulting non-convex "effective potential" still remains to be 
understood. 



In [18 1 it was pointed out that the usual Landau-gauge effective potential and the 



gauge invariant effective potential lead to different predictions for observables like latent 
4 For a discussion and references, sec |G2|. 
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heat, critical temperature etc. Our derivation of the gauge invariant effective potential 
in this section demonstrates that these differences are a consequence of different choices 
of the resummation procedure. Based on the arguments given above, the asymmetric 
resummation which treats symmetric phase and Higgs phase differently, appears better 
justified. 

The results of [IB] were obtained by performing an expansion around the tree-level 
minima in the Higgs phase and in the symmetric phase. This approach has been consid- 



ered at two-loop level |25| , P7j , and some problems have been discussed by Laine |[27|| . In 
this section we have reproduced the results of [|l^] by expanding around local minima of 
the effective potential V which includes quantum corrections. As we shall demonstrate 
in the following section, this procedure can be extended to two loops. At this level the 
fundamental infrared problem of the symmetric phase will also become apparent, and we 
shall discuss under which conditions the loop expansion can yield a good approximation 
of thermodynamic quantities relevant for the electroweak phase transition. 



5 Two-loop results 

In this section we shall extend the asymmetric resummation to two loops in order to 
examine the convergence of the perturbative expansion. In recent years several two-loop 
calculations have already been carried out. In [|J all two- loop contributions involving only 
the gauge coupling were evaluated, which yield the effective potential V to order g 4 , A. 
In this calculation scalar masses were set equal to zero. A complete two-loop calculation 
of V to order g 4 , A 2 , including scalar loops, has been carried out in @, where also the full 
standard model has been considered. Here the masses (28),(|29), with J = 0, were used 



in the symmetric phase and the Higgs phase. The same result has also been obtained 
by using as intermediate step the effective three-dimensional theory, which is obtained 
by integrating out modes with non-zero Matsubara frequencies |28j. Furthermore, the 



effective potential of the three-dimensional theory has been evaluated in general covariant 
29] and 't Hooft background gauges |30] . 



In all these calculations the resummation of scalar masses has been performed in 
the same way in the symmetric and the Higgs phase. However, on physical grounds, as 
explained in the previous section, an asymmetric treatment of the two phases appears 
more appropriate. In the Higgs phase scalar masses are given by eq. PH ) and the counter 
term is defined by the sum 5Sp + 5Sp, as discussed in sect. 4. In the symmetric phase 



If 



the scalar masses are self-consistently determined from eq. (|33] 



ml = m l = -^-V(T } 0;J) , (46) 

<p Off 

where the potential V is calculated with the counter term 8S% given in eq. (|37]) . Note, 
that in the two-loop calculation the counter term to be inserted in the one-loop graph is 
0(g 3 ), whereas the scalar mass determined from eq. (|3^) is of higher order in g. 

The two-loop potential V for arbitrary masses m a and m n can be extracted from 
[[I]. [/J. We have listed the individual contributions corresponding to the graphs of fig. 9 in 
the appendix, omitting terms independent of (p and terms which cancel in the sum. The 
two-loop potential contains terms linear in m a and m n . In the symmetric phase, where 
the counter term is given by (|37|) , these terms cancel in the sum. In the Higgs phase they 
contribute to the potential. 

In the symmetric phase the scalar masses are determined self-consistently by eq. 
fl3~3|) . With m T = g<p/2, the two-loop potential yields a contribution which diverges 
logarithmically at ip> « (cf. fl56p), 



ml*-^£- 2 T*\nf3m T . (47) 

Following H we regularize this divergence by means of a "magnetic mass" term. In eq. 
(fT7|) we substitute m\ = g 2 p> 2 /4 + 7 2 g A T 2 / (9-7T 2 ) . In the following we shall use 7 = 1, 
which follows from one- loop gap equations We have checked that the results of 

our numerical analysis change only insignificantly if we vary the parameter 7 between 
0.3 and 3.0. The deviation of the most sensitive quantity AQ/T C 4 from the plotted 7 = 1 
result is invisible for small Higgs mass and increases up to 8% at tuh = 70 GeV. Eq. 
( |33|) for the scalar masses can be solved iteratively. In the first step one inserts in the 
two-loop potential V the one-loop scalar masses (|28|) . The mass m a obtained in this first 
step of iteration is already a very good approximation to the exact solution in the range 
of Higgs masses which we shall consider. 

Given the two-loop potential in the symmetric phase and the Higgs phase, we can 
numerically determine the critical temperature T c , where the two potentials at their 
respective local minima are degenerate. Differentiation with respect to temperature and 
the external source J then yields latent heat and jump in the order parameter p = 2$t$. 
The results are shown in figs. 4-6, assuming standard model values m w = 80.22 GeV 
and v = 251.78 GeV at zero temperature. 

In fig. 4 the critical temperature in units of the Higgs mass is plotted as function of 
the Higgs mass in the range 30 GeV< mn <80 GeV. Below 30 GeV the high-temperature 
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expansion is unreliable, and above 80 GeV the convergence of the loop expansion dete- 
riorates rapidly. The two-loop results are compared with one-loop results for the "old 
resummation" and for the "new resummation" described above. In the "new resum- 
mation" the one-loop result is lowered whereas the two-loop result is increased, improving 
the convergence of the loop expansion considerably. The relative change of an observable 
from one-loop to two-loop may be characterized by 5 = 2\Oi — 2 \/(Oi + 2 ). In the 
"new resummation" 5 ~ 0.04 in the whole range of Higgs masses considered. 

The jump in the order parameter at the critical temperature T c is shown in fig. 5. In 
the case of "old resummation" (p c corresponds to the position of the second minimum, 
for "new resummation" it is \^Ap. Again the "new resummation" procedure improves 
the convergence significantly. The relative error increases from 5 ~ 0.01 at tuh = 40 GeV 
to 5 ~ 0.2 at m H = 70 GeV. 

In fig. 6 the latent heat AQ in units of the critical temperature is plotted as function 
of the Higgs mass. This is a measure of the strength of the first-order phase transition. 
Like the jump in the order parameter it decreases with increasing Higgs mass. Based on 
the Clausius-Clapeyron equation, which is satisfied exactly, we expect AQ ~ Ap. Hence, 
the convergence should be worse than for the order parameter. This is indeed the case. 
The relative error increases from S ~ 0.1 at nin = 40 GeV to 5 ~ 0.4 at m# = 70 GeV. 

It appears satisfactory that the new, asymmetric resummation procedure leads to an 
improved convergence of the perturbative expansion. The "old resummation" is based on 
a systematic expansion of the free energy in powers of the couplings g and A. Therefore 
in the symmetric phase, where the scalar masses are positive, the differences can only 
correspond to contributions of higher order. This ambiguity in the resummation proce- 
dure is analogous to the well known dependence of results in perturbative QCD on the 
renormalization scheme. However, setting the Goldstone mass to zero in the Higgs phase 
does not correspond to a simple rearrangement of the series, since terms non-analytic 
in m n are present. The correct treatment of the Higgs phase is the main qualitative 
advantage of the new resummation. 

6 Comparison with lattice simulations 

In the previous sections we have presented a quantitative description of the electroweak 
phase transition in the perturbative approach. Using different resummations in the Higgs 
phase and in the symmetric phase, we concluded that for Higgs boson masses below 
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~ 70 GeV the perturbative approach converges, and that this description is therefore 
self-consistent. However, for Higgs masses above the present experimental lower bound, 
tuh > 63 GeV, the loop expansion becomes unreliable. In this relevant range of the 
parameter space the electroweak phase transition can only be understood by means of 
non-perturbative methods. Lattice Monte Carlo simulations provide a well defined and 
systematic approach to study this problem. For all Higgs masses, non-perturbative effects 
may be important in the symmetric phase. By comparing data from lattice simulations 
for the SU(2) Higgs model with the perturbative results one can hope to identify non- 
perturbative features and to achieve a better understanding of the electroweak phase 
transition. Therefore, in this section we shall present a detailed comparison between data 



of the recent large scale four- dimensional lattice works [10], [TTJ and published results of 
the "old resummation" [[7]. For the considered Higgs masses they differ little from the 
"new resummation" results. 

Monte Carlo simulations for Higgs boson masses near or above the W-boson mass are 



technically difficult, thus in [[L0| [Ll|] the Higgs masses tuh ~ 18 GeV and m# ~ 49 GeV 
have been studied. Since in this parameter range the perturbative expansion converges 
rather well, one may expect agreement between perturbative and non-perturbative results 
with comparable accuracy. 

The g 3 , A 3//2 -potential of [7J] involves a high-temperature expansion up to order 
(m/T) 3 , which is unsatisfactory for m# ~ 18 GeV. Thus, we have included all one-loop 
contributions of order (m/T) 4 in our present g 3 , A 3 / 2 -potential. Note, that the numerical 
evaluation of the one-loop temperature integrals gives a result which agrees with the 
above approximation up to a few percent. 

For small Higgs boson masses the renormalization scheme dependence is non- 
negligible. Therefore, instead of the MS-scheme with Jx = T, we shall use the scheme sug- 
gested by Arnold and Espinosa Q, which includes the most important zero-temperature 
renormalization effects. In this scheme the correction to the MS-potential, used for both 
the one- and the two-loop results, reads 

where 

A 9 ^ V AX 9 9 4 (^ m w /a q\ 

Here v is the zero-temperature vacuum expectation value and mw is the W-boson mass 
at T = 0. 
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In [J 1] several observables have been determined, including renormalized masses at 
zero temperature (run, raw), critical temperature (T c ), jump in the order parameter 
(ifc), latent heat (AQ) and surface tension (a). As usual, the dimensionful quantities 
have been normalized by the proper power of the critical temperature. The simulations 
have been performed on L t = 2 and L t = 3 lattices (L t is the temporal extension of 
the finite-temperature asymmetric lattice). The L t = 3 results should be closer to the 
continuum values, and we therefore compare these data with the perturbative results. 
An exception is the surface tension for which only L t = 2 data exist. 

In the previous sections we have only evaluated T c , (p c and AQ, which follow from 
the free energy of a homogeneous phase. The surface tension is more complicated to 
calculate, since it involves the boundary between two phases. In perturbation theory 
this is related to the effective potential V in the region between the two local minima. 
For completeness, we include the surface tension in the comparison, although so far no 
satisfactory treatment has been achieved in perturbation theory for Higgs masses above 
~ 40 GeV. 

The statistical errors of these observables are normally determined by comparing 
statistically independent samples. The systematic errors can be estimated by the dif- 
ference between the L t = 2 and the L t = 3 data. As fig. 15 of |T]J suggests, the true 
systematic error for T c /rriH may be larger than this naive estimate. Thus, for T c /w,h we 
have doubled the above error. For the surface tension, where only L t = 2 data exist, the 
systematic error has been estimated as twice the statistical one. A correct comparison 
has to include errors on the parameters used in the perturbative calculation. These un- 
certainties are connected with the fact that neither the Higgs boson mass nor the gauge 
coupling has been determined exactly. Therefore, the perturbative prediction for an ob- 
servable is not one definite value but rather an interval, given by the uncertainties of m# 
and g. In this analysis only the statistical errors in the determination of the Higgs boson 
mass and the gauge coupling have been included. We have also neglected corrections due 





m H « 18 GeV 


m H m 49 GeV 


AQ from Clausius-Clapeyron eq. 


.0236(14) 


.00171(15) 


direct lattice result for AQ 


.0194(15) 


.00151(12) 



Table 1: Comparison of the latent heat in lattice units obtained by using the Clausius- 
Clapeyron equation and lattice Monte Carlo simulations. The data are from ref. J71|/ . 
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to the different renormalization conditions used for the gauge coupling g and the order 
parameter on the lattice and in the continuum. These corrections are expected to 
be of relative order g 2 . 

Let us first check the validity of the Clausius-Clapeyron equation for the lattice 
Monte Carlo simulations. The first line of table 1 contains AQ obtained by eq. ([ZH) 
from lattice results for Ap and m 2 H . The second line contains the latent heat determined 
directly by lattice simulations. As usual, the numbers in the parentheses denote the 
statistical errors. The values for m# m 49 GeV agree within one standard deviation. For 
run ~ 18 GeV, however, there is a small discrepancy. 

The comparison for <p c /T c , AQ/T^, cr/T c 3 and T c /m H is shown in fig. 7 {m H « 18 
GeV) and in fig. 8 (m# ~ 49 GeV). The dots with error bars represent the perturbative 
result at one-loop (g 3 ) and at two-loop (g 4 ) level. For each quantity the dashed lines 
show the region allowed by the statistical error, whereas the dotted lines include the 



systematic error as well. Note, that the values of g and mn/mw are those of JET]. We 



emphasize that, both in the perturbative calculation and in the lattice simulations, the 
surface tension is the most problematic quantity. 

For tuh ~ 18 GeV (cf. fig. 7) both the one-loop and the two-loop results are in 
good agreement with the lattice data. For m# ~ 49 GeV (cf. fig. 8) the two-loop results 
agree definitely better with the Monte Carlo data, except for the surface tension. The 
two plots may be interpreted in the following way. For small Higgs boson masses the 
perturbative approach is in very good shape, already the one-loop approximation gives 
a reliable result. As m# grows, the higher order contributions become more and more 
important, yet a two-loop calculation is still satisfactory for mj? ~ 49 GeV. In this range 
of parameters the non-perturbative features of the symmetric phase are not important 
enough to destroy the perturbative picture. 



7 Conclusions 

In the present paper a quantitative description of the electroweak phase transition based 
on a gauge invariantly coupled source term has been attempted. Our main result is 
that for Higgs masses below ~ 70 GeV thermodynamic observables can be evaluated 
with reasonable accuracy in perturbation theory. Above m# ~ 70 GeV the perturbative 
expansion breaks down, which is in agreement with previous estimates ||. 

The main technical achievement is the realization of different resummation proce- 
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dures in the symmetric phase and in the Higgs phase, which are expected to improve the 
convergence of perturbation theory. Such an improvement has been explicitly verified for 
the available terms up to two-loop order. Due to infrared divergencies of the non-abelian 
theory a cutoff is needed in the symmetric phase at two-loop order. We emphasize that 
due to the small cutoff dependence at not too large Higgs masses perturbation theory may 
work with some accuracy in the symmetric phase. The size of non-perturbative effects 
can only be determined by comparing perturbative results with fully non-perturbative 
lattice simulations. 

We have carried out such a comparison, based on results of recent simulations on 



large lattices ||11|| . The good quantitative agreement found for Higgs mass values ran ~ 18 
GeV and m# m 49 GeV is interpreted as evidence for the correctness of the present 
understanding of the electroweak phase transition. Non-perturbative effects present in 
the symmetric phase are neglected by perturbation theory, but they should contribute 
to the lattice results. We conclude that these effects can not be of major importance at 
small Higgs mass, since otherwise no quantitative agreement with lattice data could be 
observed. 

Applying the Clausius-Clapeyron equation to the electroweak phase transition a 
simple relation between latent heat and jump of the order parameter has been derived. 
Being in good agreement with perturbative as well as with lattice data, it strengthens 
confidence in the correctness of the treatment of the phase transition. 

The above arguments support the conclusion that our understanding of the elec- 
troweak phase transition has reached a quantitative level for Higgs masses up to ~ 70 
GeV. A strong decrease of the strength of the first-order transition with increasing Higgs 
mass is observed. However, a complete understanding of the process of symmetry restora- 
tion for large Higgs masses is still lacking. Other important questions include the descrip- 
tion of metastable and unstable states, relevant for the dynamics of the transition. 



8 Appendix 

In this appendix, for the convenience of the reader, all the contributions to the effec- 
tive potential from one and two-loop graphs (see fig. 9) are listed in MS-scheme. This 
calculation has already been performed in Q up to order g A and extended in []7|, |28| to 
the order g 4 , A 2 . The different contributions are presented according to the counterterm 
method of resummation, applied to the zero modes only PI. The scalar resummation 
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is left unspecified, since it is different in the symmetric phase and in the Higgs phase 
following the approach of section |5|. It is characterized by the values of the scalar masses 
m CT and m n , and the scalar counter term C s , present in V a . 

Tree level potential and one-loop corrections are collected in V\ while the two-loop 
contributions V a ... are labelled in correspondence with the diagrams of fig. 9. Linear 
mass terms, poles in e and terms proportional to i t (see ref. ||), which cancel systemat- 
ically in the final result, are not shown and the limit e — > has already been performed. 
This is the reason for the vanishing of V&, which would only contribute to the well known 
cancellation of linear mass terms of order g 3 . 
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Here evaluating the scalar one- and two-loop temperature integrals the constants 

ci « 5.4076 and c 2 « 3.3025 



(59) 



have been introduced following |yj and |32|, |33|, respectively. 
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Figure captions 

Fig.l Free energy at zero source term as a function of the temperature 
(m H = 70 GeV). 

Fig. 2 Free energy at the critical temperature as a function of the source J 
(m H = 70 GeV). 

Fig. 3 Gauge invariant effective potential (solid line) and its analytic continuations from 
the single phase regions into the mixed phase region (dashed line). The cross denotes the 
matching point. 

Fig. 4 Ratio of critical temperature and zero-temperature Higgs mass as a function of 
the Higgs mass. 

Fig. 5 Square root of the jump of the order parameter y? c = yj Ap in units of the critical 
temperature T c . 

Fig. 6 Higgs mass dependence of the latent heat AQ of the phase transition. 

Fig. 7 Comparison of the lattice data with the one-loop and two-loop perturbative results 
for m H a: 18 GeV. 

Fig. 8 Comparison of the lattice data with the one-loop and two-loop perturbative results 
for Tfijj ~ 49 GeV. 

Fig. 9 Two-loop diagrams contributing to the effective potential. 
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